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Abstract 

The implementation of gauge theories on a four-dimensional anisotropic lattice with two distinct 
lattice spacings is discussed, with special attention to the case where two axes are finely and 
two axes are coarsely discretized. Feynman rules for the Wilson gauge action are derived and 
the renormalizability of the theory and the recovery of the continuum limit are analyzed. The 
calculation of the gluon propagator and the restoration of Lorentz invariance in on-shell states 
is presented to one-loop order in lattice perturbation theory for SU(N C ) on both 2+2 and 3+1 
lattices. 
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I. INTRODUCTION 



The anisotropic lattice is a popular tool in numerical simulations. The usual approach 
is to make the temporal lattice spacing fine while keeping the spatial directions relatively 
coarse. The extra temporal resolution in this 3+1 discretization scheme is useful when 
measuring two- and three-point correlation functions of particles which decay rapidly. This 
is particularly useful for heavy hadrons and glueball states,for which the signal to noise 
ratio decreases rapidly with time. The glueball spectrum Jjl 0] was an early success for 
this approach and,more recently, heavy quark systems have been studied using anisotropic 



lattices 
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In this paper we consider a generalization of the anisotropic approach to include lattices 
with a 2+2 discretization. The temporal and one spatial direction are made fine, keeping 
the remaining two spatial directions coarse. The motivation for this choice is to explore the 
feasibility of calculating decays which produce high-momentum daughter particles. These 
include the phenomenologically interesting transistions B — > nlv and B — > K*j. CKM 
matrix elements are determined from such exclusive decay processes through a combination 
of experimentally measured branching fractions and theoretical calculations of form factors. 
The light daughter hadrons in these decays have non-zero momentum and so in a lattice 
calculation of the nonperturbative form factors there are cutoff effects proportional to this 
momentum in units of the lattice spacing. Therefore B — > -nlv has discretisation errors 
proportional to ap^ where a is the lattice spacing. However, the range of momenta reliably 
reached by current experimental and lattice data do not overlap. Lattice calculations work 
best with p n < lGeV but the bulk of experimental data for the exclusive decay, B — > tiIu 



lie at > 1.5GeV. To avoid model-dependence kinematic cuts can be app 



the range of lattice data to momenta where the calculation is reliable 



ied, restricting 



151 ] . awaiting im- 



proved experimental results. Alternatively, the lattice data can be extrapolated to match 



experiment 



17L Il8j ] but this introduces model-dependence and increases substantially the 



systematic error in lattice calculations and therefore in \V u b\. For B — > K*^ the situation 
is more acute since this decay happens at the maximum recoil momentum of the daughter 
meson, far from the range accessible to current lattice calculations. Therefore it is neces- 



sary to extrapolate well outside the range of reliable data 



19, 20]. once again introducing 



model-dependence and increasing the systematic errors. For this reason it has not been 
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widely studied using lattice methods and calculations to date have used isotropic lattices. 
In both cases the difficulty for lattice calculations is that errors proportional to the momen- 
tum of the order ap, grow quickly. These must be controlled to make phenomenologically 
relevant lattice calculations. The advantages then of the 2+2 lattice are twofold. Firstly, 
the fine temporal lattice spacing serves the same purpose as in the 3+1 case: correlation 
functions should be accurately determined while keeping computational costs modest. Sec- 
ondly, making one spatial direction fine and injecting all momentum along that direction 
keeps discretization errors of 0(ap) small for high momenta. 

The transition B — > ttIu has been calculated using a 3+1 discretisation scheme and the 
improved resolution in the temporal direction led to higher- momentum (0 < p n < 1.5GeV) 
particles being reliably simulated j^]. Both this calculation and the isotropic calculations 
which were reviewed in Ref. |22j find that one of the largest systematic errors in the range 
of accessible momenta is due to the chiral extrapolation. The 2+2 discretization does not 
address this issue but it is hoped that it will further extend the range of momentum available 
to lattice calculations and that statistical precision will be enhanced. 

The paper is organized as follows. In Section |H] a general anisotropic formulation for 
SU(N C ) Yang-Mills theory in four dimensions is presented. The differences between 3+1 
and 2+2 discretisations are discussed in terms of the lattice symmetries and the parameter 
tuning required. Section ITTT1 contains the framework of the perturbative calculation, already 
outlined in Ref. |3- Although the goal is a 2+2 discretization, the Feynman rules and 
the analytic procedure given are completely general and allow an exact treatment of the 
calculations. In Section IIVI we present our results and include a comparison with other 
work. The challenging algebraic manipulations were performed on a computer using a 
symbolic code which handles the dependence on the anisotropy analytically. The core of the 



code is similar to that used in Ref. 



24 



while the treatment of the lattice integrals follows 
from Ref. |2^|. Section contains our conclusions, while some technical details are given in 
Appendices El and El 

II. YANG-MILLS THEORIES ON A 2+2 ANISOTROPIC LATTICE 

In this section, formulations of SU(N C ) Yang-Mills theory on general orthogonal lattice 
types are considered. A generalisation of the Wilson action for SU(N C ) Yang-Mills theory 
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on the lattice is given by 

Sw = J>„ „ (l - ^Tr(P, u (x) + Pl(x))) , (1) 

where (3 = 2N c /g 2 with g 2 , the lattice coupling constant and N c the number of colors. P^ v 
is the plaquette in the (//, v) plane; 

P„u{x) = U„(x)U„(x + p)Ul(x + u)U}(x). (2) 

The six coefficients, in Eq. (JJJ will describe the anisotropy class. 

If some sub-sets of these parameters are given identical values, the lattice action may 
have symmetries under the interchange of axes. In particular, if the two identities 

C12 = c 23 = c 3 i(= c cc ) and c 41 = c 42 = c 43 (= c cf ) (3) 

are imposed, the lattice action is symmetric under the cubic point group, and the 3+1 
anisotropy class is realised. Similarly, if the identity 

C13 = cu = c 23 = c 24 (= C e f) (4) 

is imposed, the lattice action transforms trivially under elements of the group C 4j , <8> C± v , 
with the first constituent group comprising the rotations and reflections in the (1,2) plane 
and the second being those in the (3,4) plane. This anisotropy class is denoted 2+2. 

For phenomenologically motivated reasons, as discussed in Section HI consideration will be 
restricted to these cases where only two distinct lattice spacings are permitted. The coarse 
and fine lattice spacing will be denoted a c and a/ respectively. In this paper the overall scale 
a is chosen to be a,f. The implementation of a 3+1 anisotropy, in which the lattice spacing 
in the temporal direction is made fine keeping the spatial lattice spacing coarse, has been 



widely discussed in the literature 
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In the continuum, the corresponding sub-groups of 0(4), the rotations in Euclidean space 
that are generated by demanding the same symmetries in the coefficients of the continuum 
dimension-four operators, are 0(3) for the class described in Eq. (JHJ and 0(2) <g> 0(2) for 
Eq. (jlj). If the operators which transform trivially under these symmetries are enumerated 
there is one operator for the Euclidean group, (Tr F^F^), two operators for the 0(3) 
spatial rotation group (most conveniently denoted Tr EiEi and Tr BiE>i for chromoelectric 
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and chromomagnetic) and three operators for the 0(2) ® 0(2) case (Tr FfpFff, Tr Ff c Ff c 
and Tr F cd F cd ). 

The recovery of Euclidean invariance in the low-energy physics of an anisotropic lattice 
requires a parameter tuning, in contrast to the isotropic case. Additional differences arise 
between the 3+1 lattices, for which a one-parameter tuning suffices and 2+2 lattices which 
necessitate a two-parameter tuning. The importance of parameter tuning, in particular for 
the 2+2 case is emphasised in subsection III Al 

Finally, it is interesting to note that taking the anisotropy £ — > oo yields, up to a gauge 
transformation, the Hamiltonian limit of the theory on a 3+1 lattice. For the 2+2 anisotropy 
class, this limit can not be taken naively; a theory in which only two of the four dimensions 
are discretised is not regularised. 

A. 2+2: the need for tuning 

As the continuum limit is approached for the 3+1 lattice theory, the physics of the 
Euclidean invariant Yang-Mills theory is reproduced, provided care is taken to account for 
the different grid spacings, a c and a. This can be achieved by tuning the relative weights 
of the two sets of coefficients in Eq. Q to ensure the ratio of scales, £ = a c /a measured 
by a physical probe takes its desired value or alternatively, an arbitrary choice of the two 
parameters c cc and c c f can be taken and the ratio of scales measured post-hoc. The 3+1 
anisotropic theory is certain to be in the same universality class as the desired theory since 
the single free parameter (the relative weight of c cc and c c f) determines the ratio of scales £. 

For the 2+2 case, an important distinction arises; while there are only two distinct 
dimension-four continuum operators in the 3+1 class, there are three for the 2+2 case. As 
a result, the recovery of a Euclidean invariant continuum theory is not guaranteed since 
there are two free parameters and only a single ratio of scales. As a result, the general 2+2 
lattice theories can lie in a larger universality class than the continuum four- dimensional 
Yang-Mills theory. 

For this reason, care must be taken to ensure the recovery of Euclidean invariance in 
the continuum limit for a 2+2 simulation. The relative weights of the three operators in 
the action must be determined to ensure Lorentz invariance in on-shell Green's functions. 
This tuning can be achieved perturbatively or by restoring symmetries in a non-perturbative 



5 



calculation of, for example the static inter-quark potential. In the remainder of this paper, 
perturbation theory is used to determine the parameters in the action. A paper describing 
the nonperturbative tuning of the parameters is in preparation j^j. 



B. c^ u from perturbation theory 

At tree level, it is straightforward to determine the values of the coefficients in the action. 
For the 3+1 lattice action, they are c cc = l/£, c c / = £ and for the 2+2 case, they are 
c cc = V£ 2 , c cf = 1 an d Cff = £ 2 . The one- loop definitions of c^ v can then be parameterised 

(jl, v fine 

fi coarse(fine), v fine(coarse) (5) 
/i, v coarse 

a coarse(fine), v fine(coarse) 

(6) 

/i, v coarse 

In both cases, an overall multiplicative weight can be absorbed into a redefinition of the 
field integration variables, and this is redundant in the action-tuning procedure since a 
symmetry is being enforced in on-shell Green's functions. This change of variables will 
become important once matrix-element matching of gluon fields is being performed. 



as 



and 



e (l + 4 ] g 2 + 0(g*)) 

^ (2+2) = J i (i+^y+oG? 4 )) 

£ (l + r£V + OQ? 4 )) 

,(3+1) _ j * a+^SV+ote 4 )) 



III. PERTURBATION THEORY 



We can now proceed to discuss Feynman rules for a general four-dimensional anisotropic 
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Wilson-like action, already outlined in Ref. |23j. We follow the notation given in Ref. 
and the interpretation of the anisotropy as a difference in momentum cut-offs will be our 
guideline. The matching of the lattice gluon action in Eq. (0) with its continuum counterpart 
is made clearer by expressing the link variables U^n) = e l ^ n ' in terms of dimensionless 
fields 

Un)=^{n)T\ (7) 

where T b are the SU(N C ) generators in the fundamental representation satisfying the rela- 
tions, [T a ,T b ] =if abc T c and Tr (T a T b ) = ^5 ab . The dimensionful gluon fields, A^x) can 
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then be reintroduced in Eq. (JJJ) using the relation 

^(n) = ga^Al(x), (8) 

where 

{1 n fine , x 

(9) 
£ \i coarse 

is the anisotropy index. By taking into account the Jacobian 

n x 

where d denotes the number of coarse directions, one easily gets Eqs. (0) and (jHJ). The 
continuum limit of the Eq. (JIJ is then written as 

Sw = Y, \c^F> v {n)Fl{n) + O = 1 £ F>„(x)F» u (x) + 0(a 2 ) , (11) 

where (9 is an irrelevant operator in the continuum involving terms with three or more gluon 
fields from which the non-trivial contribution to the Feynman rules will arise, while 

FlM = d^Hn) - d^%n) - gf hc %{n)<j> d v {n) , (12) 

and 

d?<p(n) = <p(n + p,)-<p(n) (13) 
are the dimensionless field strength and lattice right derivative, whereas 

= d >^ x ) ~ d * A » ~ 9f bcd Al(x)At(x), (14) 

and 

d*<f{x) = — {<p(x + aCnP)- <p{x)) (15) 

are their dimensionful equivalents. 

As is well known, the Wilson action has 0(a 2 ) discretization errors in the evaluation 
of various physical quantities. In addition, the lattice regularization gives rise to finite 
renormalization coefficients when compared to other continuum schemes. These effects can 
be reduced by adding irrelevant operators in the action which reduce discretization artefacts, 
or by improving the convergence of renormalization coefficients to better match continuum 



quantities 
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3 it ] . In any of these cases an exact and completely algebraic treatement of the 



Feynman rules is always viable, even at orders higher than one As for the 

lattice integrals, while a completely numerical evaluation presents no problem at first order, 
any one-plaquette action will change the structure of the 0(a 2 ) (irrelevant) terms in the 
four and higher gluon vertices but will not change the propagators, the three-gluon vertex, 
the measure or the gauge-fixing terms, which are all fixed by the Haar measure and the 
naive matching with the continuum limit 3]| • This means that once the first-order analytic 
technique for the Wilson action is developed it can be applied to any other one-plaquette 
action. For gaug e actions with Wilson loops which extend two grid spacings, the technique 
given in Ref . |2ja (for fermion propagators) can be adapted while for actions with even larger 
loops, suitable techniques can be developed in the same spirit. Using coordinate-space 



38l |39|. Finally, mean- 



methods, analytic results are also available for higher orders 
link improvement schemes |4(| are straightforward to implement once the corresponding 
Feynman diagram is written down. 



A. The gauge fixing and the gluon propagator 

For perturbation theory it is necessary to choose a gauge fixing and proceed to define the 
other parts of the action needed to write the Feynman rules. In analogy with the isotropic 
case we choose the gauge-fixing term such that 



(16) 



where dj^ip(x) = l/(a^) (<f(x) — ip(x — a£ M /t)). It is straightforward to show that 

s gf = 4 Tr ^ 2 for T = , (17) 



where now d^ip(n) = <p(n) — ip(n — fi) and 



e- 1 



(j, fine , 
// coarse , 



(18) 



gives the desired form once a = £~ 2l+d ~ 2 a. The choice I = (d — 2)/2 would be easiest but 
we shall see in the following it is not the most convenient. The gluon propagator can now 
be calculated from the two field terms in Eq. ()12|) and the definition of the gauge fixing in 
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Eq. (|17|). The functional form is similar to the isotropic case, namely 

k 5 ab f k k \ 

^ a u > b = |J ( s n» ~ 0- ~ a)-^ I > ( 19 ) 

where 

L = 4" sm(^U , P = V P , k „ e 4" V (20) 

are the dimensionful lattice momenta spanning the anisotropic Brillouin zones. All the 
anisotropy is now encapsulated in this form of the Brillouin zone. This will prove to be a 
constant pattern: for any vertex or propagator which has a continuum analogue, the form 
in terms of the (anisotropic) dimensionful momenta will resemble the isotropic ones, since 
the asymmetry now lies solely in the different momenta cut-offs. 



B. The Haar measure 

On the lattice the functional integral is obtained by integrating over the dimensionless 
gauge links, but to do perturbation theory it must be expressed in terms of the dimensionful 
gauge fields. The Jacobian resulting from the expression of the Haar measure in terms of 



Eq. © is 



J — \ det 



-MMt , M(0») = , 0»=«^(n)t<\ (21) 

1 J rfM 



where the 4>^{n) are the same as in Eq. (J7J) and the t a are the adjoint generators of SU(N C ). 
Reexpressing J as e _5mcas gives 

2(1— cos (bJn) 



mens 



1 i z i ± - cut, I \ 

^-V 1 !' (22) 



When expanded to lowest order in the fields (which is all that is needed for one-loop two- 
point function calculations) this reads 

S^*^^N^Mn)fi(n), (23) 

which putting back the dimensionful fields gives the vertex 

k k' a 2 N £ 2 

W ^ b = -(2n) 4 6\k + k')^6,J ab ^. (24) 

By convention all gluon momenta are incoming. 



C. Fadeev- Popov ghost fields 



The Faddeev- Popov determinant, which forces the integration only on a section of the 
gauge orbits, must also be included in the action. Using T in Eq. (jT7jl to enforce the gauge 
condition on M and using the response of the gauge fields on the lattice to (infinitesimal) 
gauge transformations [29] 

DM = M-\$ ll )d2 + i$ li , (25) 
where and the inverse adjoint- valued operator, M are given in Eq. ([21)1 . we obtain 

3* = - E ?{n)tf>%Bt(<l>)*{n) , (26) 

n,fi 

where c and c are the dimensionless lattice ghost fields, introduced to make Eq. local. 
From the expansion of M _1 

M- 1 (^)=U + t -$»-^l+... (27) 



x 



and by reintroducing dimensions, Eq. (|26|) can be written as 

S gh ~ -^-^^(a;)rt^y(x)+^ c ^^)^[^(x)(l + ^^)]c 6 

+ i^^^O^* d } a6 (^c a (x))(^c 6 (x))^(x)^(x; 

by including only terms which are relevant for one-loop two-point functions. The factor 
at the front of S^h gives rise to spurious coefficients with no continuum analogue which 
however, cancel in any graph with no ghost outer legs. Alternatively the coefficients can 
be absorbed in a redefinition of the ghost fields, £^ _1 " 2 c — > c or by setting I = 1, thus 
rescaling the lattice gauge-fixing parameter a = £ 4_d a. For each of these alternatives the 
ghost propagator and the two ghost one- and two-gluon vertices are 
k „„, 1 



a ■ 

p',b 



■ b = 5 ab — (28) 
k 2 

V ^kl^' c = ig(27r)^(k+p-p , )f ab %p l ,, (29) 



4 

p, b p', a 

\ 4 1_ 

12 



g 2 a 2 (2n)*5\k + k' + p - p'){t c , t d } ab 5, u ftp^, (30) 



k,/j,,c k',p,d 
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with = cos \p^a^^. As usual the vertex which has a continuum analogue carries no 
explicit dependence on the anisotropy. 



D. Gluon vertices 



Taking into account the three-gluon terms arising from Eq. |Q, one of which comes from 
the irrelevant operator O defined in Eq. flllj) gives 

S§? = i A » + ^My*)) (d*Al(x))Al(x), (31) 

X ^ ' 

which, after a Fourier transformation, leads to an expression equivalent to the isotropic case 
k, n, a 

"^5^ k", A, c = ig{2ixY5\k + k' + k")f abc (s uX (]f^k') ~k v 
k '> "> b + 5^{k^k") v k\ + 5^(k^k) x k%) . (32) 



Taking the four-gluon terms arising from Eq. ([T|. mostly coming from (9, (see 

Ref.t2Sl 

for 

details), and rechecking the cancellations which must occur due to Bose symmetry, one finds 
the four-gluon vertex whose expression is given in Appendix lAl 

E. One-loop correction vertex 

The one-loop corrections to the action coming from Eqs. (jHJ) and (JHJ) give rise to an extra 
vertex which reads 



where 



„(i) 



r)ff /i, i/ fine 

77^ /i coarse(fine) and z/ fine(coarse) (34) 



(i) 

rfcc /i, f coarse 



In the 3+1 case one can either set rfff to zero or leave it free, as by the traceless property 
of the vertex it will always cancel. 
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F. The continuum limit and anisotropic renormalization 



The calculation of the one-loop corrections to the gluon self-energy, E^a^p), using the 
Wilson action involves five Feynman diagrams j^j and Eq. (|33|) . Each diagram is a function 
of the external momenta p and can be written as 

d A k 



G(p) 



(2*) 



-F(k,p) 



(35) 



where k is the integration momenta. Since we are interested in the continuum limit of 
Eq. (JHHJ), if the integral is ultraviolet convergent we can simply substitute the function 
F(k,p) with its continuum equivalent. Otherwise, if Eq. ()35|) is divergent and contains only 
massive propagators so that F(k,p) is finite for any set of momenta going to zero, one can 
use the lattice version of the BPHZ technique 0, 

r v J* 1 

G(p) = 



d 4 k r x 1 



44j | by writing 
d d 



n=0 



dp, 



., F(k,p) 
Vp,,,, ■' i- -')■ 



(36) 



f d 4 k ^ 1 ( d d 

V (27r) 4 ^n! V9p Ml % n 



P) 



p=0 



G c (p) + G L (p) 



where np is the degree of divergence of the diagram. The first term in Eq. is ultraviolet 
finite j3, |3| and therefore its continuum limit can be taken, making it independent of 
£. All the effects of the lattice regularization remain in the second term, which is simply 
a polynomial in the external momenta with coefficients given by zero-momentum lattice 
integrals. 

On the other hand, if a diagram contains massless propagators, as in our case, more care 
is needed: indeed an expansion around p = can give rise to infrared divergences. A simple 
recipe is to introduce an intermediate infrared regularization. Given an anisotropic cut-off 
the introduction of a mass, m in the propagators is the most suitable solution |25j, 3.2, 45. |, 
while dimensional regularization is popular in the literature for the isotropic case 3J3, l41 |. 
G c (p) and G L (p) are then divergent for m — > separately but the divergences cancel in the 



sum. 
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IV. RESULTS 



In this section we present the results of the one-loop correction to the gluon propagator 
using the Wilson action in the Feynman gauge for a general anisotropic lattice in four 
dimensions. Applying the procedure explained in Section (|1V Aj) we calculate the values of 
the coefficients which restore on-shell Lorentz invariance, where the definition and treatment 
of the lattice integrals £>g is given in Appendix JBJ). 



A. On-shell Lorentz invariance 

As explained in Section m the one-loop propagator obtained from the Feynman rules 
given in Section IIHI will not in general satisfy Lorentz invariance. The free parameters in 
Eq. (}34|) must be tuned to restore the symmetry j3, 22, 4^. From Eq. (}33j) it is clear that 
not every Lorentz-breaking term can be cancelled, as expressions of the form 8^ /a 2 and 
dfivPfi arise from each diagram. The first non-trivial result of our calculation is that, for any 
value of d and £, such terms cancel exactly when the diagrams are summed up, just as in 
the isotropic case, independently of the tuning procedure. 

The remaining Lorentz-breaking artefacts arising from the mixing of longitudinal and 
transverse field components have the correct structure and can be cured by tuning. We 
choose to fix the parameters in the action by demanding the recovery of Lorentz invariance 
for on-shell physical soft gluon modes, imposing that the two physical eigenvalues of the 
one-loop propagator vanish for E 2 = p 2 . Since this is a gauge-invariant condition we restrict 
ourselves to a particular gauge. We have chosen the Feynman gauge for which a = 1. By 
injecting the gluon momentum in any possible direction and calculating the eigenvalues and 
eigenvectors of the propagator we obtain a general condition for rj^ — r$ and r]f) — ry^f , 
independent of the direction of p^. The value of remains unconstrained as it can be 
reabsorbed in a one- loop /3-shift. 

To identify the physical eigenmodes, the momenta directions with residual symmetry 
are first identified; In the case of the 2+2 lattice, this corresponds to any momenta in the 
coarse-coarse or fine-fine planes. Then for more general momenta, the axis of propagation 
was smoothly varied away from these symmetric cases, and the eigenmodes continuously 
traced. This investigation lead to a generalisation of the polarisation condition, p^e^ = 
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for on-shell gluon polarisation vector e„. The resulting 2+2 lattice polarisation condition is 



PuZ^u^u — 0, 



(37) 



with Z^v is the diagonal matrix = 5^ v (Z c 5^ c + ZfS^f) and Z c , Zf the coarse and fine 
gluon field renormalisation coefficients. S^f and <5 MC are one(zero) if fi is fine(coarse). 

B. The one-loop coefficients for the restoration of Lorentz invariance 



For 2+2 anisotropic lattice the one-loop coefficients for the restoration of Lorentz invari- 
ance are 

1 



„(i) _ „« 

'Ice 'lef 



2N r 



£|(1,1) 



+ N C 



16 6 V2 



1 

e 



B|(i,i) 



BU2, 1) 



5 11 \ 61(2,1,1)^1 
2 -I I ) -I i- 

2£ 4 2£ 2 / 6 



(38) 



and 



(i) (i) 

Iff ~ Vcf 



2N r 



1 1 

4 ~ 2tf 



+ 



1/1 1 \ „ / 1 

l(2-e) +W (2 + 3^ 



1) 



^(2,l,l)i 



2 ^ 4 ^ ' 'V3£ 4 £ 2 6/ 3£ 2 

which diverge logarithmically with £, while for 3+1 we have 



(39) 



„(i) _ J 1 ) 

'Ice 'lef 



Nc 



B C J1,V. 



B € (l) (i 2 19 7 \ ~g 
6 V 3 6 2£ 2 / 4 



1 



1 



^ 1 )(^^ + ^)-l]-2^[^M)(i + ^)-i 



(40) 



46. The 



which agrees with the general result of Ref. bjj and the N c = 3 result of Ref. 
techniques used by both groups involve extrapolations in some suitable parameter, absent 
in our treatment. The functions B^ are defined in Appendix El Using the results given 



there t 



re Hamiltonian limit, £ —>■ oo can be treated analytically also and agrees again with 



Ref.E. The difference in the parameters, rjcV — Vtf an d vff ~ Vtf appears as a polynomial 
in N c with terms proportional to 1/N C and N c only, as expected. Numerical values for the 
coefficients in this polynomial for a range of anisotropies are given to ten decimal places in 
Tables (HO. 
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C. The gluon self-energy in the Feynman gauge 



Once the values in Eqs. (JHHJ), an d PHj) are inserted, the one-loop correction to the 
gluon self-energy reads 



^(i) 



For 2+2 we have 



.4 



(p) = g 2 5 ab (5,„p 2 - p^) {A - t$ - B [S^S VD - <W<W) } • 



(41) 



8iV, 



AC- 



1 gg(l)/l 1 x gg(2,l) / 7 
16 2 V3 4£ 2 / 4£ 2 V3£ 2 



^(2, 1, 1) 
12 



♦fx® 



B = N r 



while for 3+1 



12 



+ 



2£ 2 



4e 2 v ^e 2 



(42) 
(43) 



2N£ 



g|(M) l 

e 2 



1 | 5^ ^(1,1) *%(2,l) /7 1 

8 3 4£ 2 4£ 2 V£ 2 3 



&(D/i 



6 V3 2^ 



5 



gg(l)/l lx gg(2,l) / 3x 

6 V2 e 2 / 4^2 v ^ ^; 



where 



J^(-io g (aV)+F (e)-7B+f : 

l07T z 9 



(44) 
(45) 

(46) 



We stress that the term proportional to B in Eq. (|41|) does not break Lorentz invariance 
but arises from the difference in the renormalizations of the fine and coarse fields A^. This 
difference must be taken into account when calculating physical quantities such as the A 
parameter or matrix elements. After including the difference, the physical modes will still 
be transverse. 

Setting £ = 1 we recover from Eq. (J41j) the well-known value of the gluon self-energy 

mm 



3W(p) = 5 ab (W 2 -^)[l-^{ T ^(-^og(pV) + 



28 
~9~ 



IN, 



A;( 72^48^ 



5 ,n N 1 

(Fo-7^) + ^ 



(47) 



where Z and F are standard integrals in perturbation theory on the isotropic lattice, and 



are defined in Ref. 



45l In order to calculate coefficients at the one-loop level in mean-link 
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improved perturbation theory, the evaluation of either the plaquette, (TrP^) or the link 
trace in Landau gauge, (Tr U^) is required. We find 



and 



N, 



1 



2^ 



(48) 
(49) 
(50) 



N, 



(TrU c ) 



1 



^ e (l)-B|(2,l) 1 



B € (l)-^(2,1) (l-£) 



(51) 
(52) 

where B^(q, 1) satisfies Eq. ()B4|) . Numerical values of the iV c -polynomial coefficients, rjcl\ 
rjty and ryj) for a range of anisotropies and including these improvement terms are given in 
Tables (jl I III VIVIVIjl . 





(i) (i) 

Vcc - %f 








l/N c N c 


l/N c 


N c 


1 











2 


-0.1161042701 0.0958601237 


0.0775690169 


-0.0741548219 


3 


-0.1886258210 0.1529850132 


0.1010941459 


-0.1044390683 


4 


-0.2391230676 0.1920838712 


0.1107973557 


-0.1217185361 


5 


-0.2774799637 0.2215115752 


0.1156439679 


-0.1335092263 


6 


-0.3083219799 0.2450412380 


0.1183899089 


-0.1424129673 


7 


-0.3340937184 0.2646297423 


0.1200891686 


-0.1495669092 


8 


-0.3562244833 0.2814071688 


0.1212112364 


-0.1555537098 


9 


-0.3756169823 0.2960809268 


0.1219898822 


-0.1607076435 


10 


-0.3928761275 0.3091218866 


0.1225517876 


-0.1652371442 



TABLE I: A^-polynomial coefficients of 77, 
anisotropic lattice. 



(i) 

cc 



7 (1) 



and 



(i) 



'Iff 



rfcf for £ 



10 on a 2+2 
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l/Nc N c 

1 

2 -0.0853420430 0.0663515978 

3 -0.1202379609 0.0892976526 

4 -0.1389282597 0.1007025039 

5 -0.1504977336 0.1074686585 

6 -0.1583451793 0.1119335310 

7 -0.1640119208 0.1150960459 

8 -0.1682939409 0.1174518583 

9 -0.1716426764 0.1192740369 

10 -0.1743328788 0.1207251404 

TABLE II: 7V c -polynomial coefficients of r^c — -qf) for £ = 1, ■ ■ ■ , 10 on a 3+1 anisotropic lattice. 





(i) (i) 




_ J 1 ) 




l/Nc N c 


l/Nc 


N c 


i 


-1/16 1/16 


1/16 


-1/16 


2 


-0.2366564052 0.2164122588 


0.1981211519 


-0.1947069547 


3 


-0.3454387314 0.3097979237 


0.2579070564 


-0.2612519788 


4 


-0.4211846015 0.3741454050 


0.2928588896 


-0.3037800700 


5 


-0.4787199456 0.4227515571 


0.3168839498 


-0.3347492082 


6 


-0.5249829699 0.4617022280 


0.3350508989 


-0.3590739573 


7 


-0.5636405777 0.4941766016 


0.3496360278 


-0.3791137684 


8 


-0.5968367250 0.5220194105 


0.3618234781 


-0.3961659514 


9 


-0.6259254734 0.5463894180 


0.3722983733 


-0.4110161346 


10 


-0.6518141912 0.5680599503 


0.3814898513 


-0.4241752080 



TABLE III: ^-polynomial coefficients of - r]^ and rf$ - r$ for £ = 1, • • • , 10 on a 2+2 
anisotropic lattice using mean (coarse) plaquette improvement. 
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(i) (i) 

'CJ 


„W. 
'Iff 
ii 


_„(!) 
c j 




1/N C N c 


1/N C 


N c 


l 


-3Z /8 3Z /8 


3Z /8 


-3Z /8 


2 


-0.2304051667 0.2101610203 


0.1918699134 


-0.1884557162 


3 


-0.3391678963 0.3035270885 


0.2516362213 


-0.2549811436 


4 


-0.4150758989 0.3680367024 


0.2867501870 


-0.2976713674 


5 


-0.4727526305 0.4167842420 


0.3109166347 


-0.3287818930 


6 


-0.5191216331 0.4558408912 


0.3291895621 


-0.3532126205 


7 


-0.5578576084 0.4883936322 


0.3438530585 


-0.3733307991 


8 


-0.5911125215 0.5162952070 


0.3560992746 


-0.3904417479 


9 


-0.6202462174 0.5407101620 


0.3666191173 


-0.4053368786 


10 


-0.6461700016 0.5624157607 


0.3758456617 


-0.4185310183 



TABLE IV: iV c -polynomial coefficients of t] cc — rj^J and ryJ — rf^J for £ = 1, • • • , 10 on a 2+2 
anisotropic lattice using mean (coarse) link improvement. Zq is defined in Section [IV Bl 



l/Nc N c 

1 -1/16 1/16 

2 -0.1694788603 0.1504884150 

3 -0.2118450434 0.1809047350 

4 -0.2337180290 0.1954922732 

5 -0.2468908748 0.2038617997 

6 -0.2556482395 0.2092365912 

7 -0.2618777619 0.2129618870 

8 -0.2665309579 0.2156888753 

9 -0.2701369251 0.2177682856 

10 -0.2730125217 0.2194047833 

TABLE V: A^-polynomial coefficients of 7$ — rf£) for £ = 1, • • • , 10 on a 3+1 anisotropic lattice 
using mean (coarse) plaquette improvement. 
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£ 


„« _ 

'Ice 






l/N c 




l 

2 


-3Z /8 
-0.1643550433 


3Z /8 
0.1453645980 


3 


-0.2068380749 


0.1758977666 


4 


-0.2288323401 


0.1906065843 


5 


-0.2420838939 


0.1990548187 


6 


-0.2508913882 


0.2044797400 


7 


-0.2571539394 


0.2082380645 


8 


-0.2618297768 


0.2109876942 


9 


-0.2654518384 


0.2130831989 


10 


-0.2683392407 


0.2147315023 



TABLE VI: iV c -polynomial coefficients of rj^ — r$ for £ = 1, • • ■ , 10 on a 3+1 anisotropic lattice 
using mean (coarse) link improvement. Zq is defined in Section fl V Bl 
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V. CONCLUSIONS 



In this paper, a generalization of the Wilson discretization to an anisotropic lattice with 
two coarse and two fine directions has been described. In particular, an important distinction 
between the 2+2 anisotropy and the well-established 3+1 case has been emphasized. The 
difference is that the coefficients in the 2+2 action must be determined before simulation to 
ensure Lorentz invariance, while any 3+1 action leads to a Lorentz-invariant theory once the 
ratio of scales, £ is determined. We are currently investigating non-perturbative techniques 
for computing these coefficients in the 2+2 case 

The main result of the paper was to compute these coefficients to first order in pertur- 
bation theory. While the focus of the calculation was on determining the Feynman rules for 
the 2+2 lattice Wilson gauge action, a more general prescription was developed to allow the 
3+1 case to be investigated as well. This allowed us to check our results against previously 
published work. The results in a mean-link improvement scheme were presented. 

The usefulness of this scheme arises from the need to make accurate calculations of 
form-factors and matrix elements at high momentum. This paper establishes the tools for 
perturbation theory calculations on the 2+2 anisotropic lattice, which will be important later 
when computing the matching factors to link calculations of weak-decay matrix elements 
to their continuum counterparts. As part of this program, quark fields on 2+2 anisotropic 
lattices are under consideration. 



APPENDIX A: THE FOUR GLUON VERTEX 




/i, a 



s,p,d 



-g 2 E f^f^ 



a 



pp^vX 
2 



(q - r)^k - s) u - —i^Kq^r^ 
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+ (b <-> c, v «-> A, q <-> r) + (b <-> d, v <-> p, g <-> s)} 



e 



APPENDIX B: BOSONIC INTEGRAL EVALUATION 

Apart from the denominators in Eq. (|19j) and (|28|1. all the £ dependence in the evaluation 
of Feynman diagrams is polynomial in £ and Having changed integration variables to 
the asymmetric Brillouin zone, only lattice zero-momentum integrals of the form 

B,(q, n u n 2 , n 3 , n 4 ) = ^ (27r)4 Dfl(fc>m)g ( B1 ) 

remain. In Eq. (jBl}, g and rij are positive integers, fc M = 2sin(A; M /2) and the £ dependence 
is in the inverse bosonic propagator 

D B (k,m) = J2K + l2^K + m2 - (B2) 
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45. In the 



alone. To evaluate these integrals, we adapt a technique given in Refs. 
following, when one of the rij is zero it is omitted as an argument of B^, while, where no 
confusion can arise, the index in is dropped and £>|, <6| or Bp c denotes integrals whose 
numerators have a fine, coarse or mixed momentum. In Refs. |25j, |45j, a set of recursion 
relations was defined to reduce every relevant integral to a linear combination of three basic 
integrals. Although similar relations exist in the asymmetric case, the lack of a complete 
symmetry among the indices, n« makes this reduction more difficult. The relations can still 
be used to reduce the number of integrals which must be calculated considerably and to 
prove the cancellations of Lorentz-breaking terms. A first set can be obtained by expanding 
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the trivial identity 

B 5 (q -n,n a ,...,n b )= / 

J —n 

which gives relations of the type 



* d 4 k k 2 a n «---k 2 b nb D B (k } m) n 



D B {k,m)<> 



d 



(4 - d)B((q, 1) + ^Bl(q, 1) = Bt(q - 1) - m 2 ^(g), 



(B3) 



(B4) 



d{ - d Z l) B\{q, 1, 1) = B^q - 2) - 2m 2 B^q - 1) + m%(q) 



(4 - d)B((q, 2) + (4 - d)(3 - d)B((q, 1, 1) + ""^ "' iff 8 (g, 1, 1) + ^(g, 2) + 



(B5) 



and so on. Furthermore, when r > 1, using the identity 



2\r 



4 W 



2\r-l 



sm /e 



D B (k } m)i D B {k,m)i ' *** q - 1 
and integrating by parts, we obtain the relation 

B c (q, . . . , r M , . . .) = J^(g - 1, . . . , r M - 1, 

+ 4B € (g,...,r Ai -l,...), 



D B (k, m 



At- 6 



is-i' 



(B6) 



— r ^(g-l,...,r,-2,...) 

(B7) 



which depends on the index fi = f,c. A final set of relations is found by using the trivial 
fact that the numerator cannot have more than four different arguments, e.g. 

Bt{q, 1, 1, 1, 1) = (d - 4)Bs(q + 1, 2 f , 1, 1, 1) + + 1, 1, 1, 1, 2 C ) + m 2 £ e (g + 1, 1, 1, 1, 1). 



:bs) 



The 2+2 case has the bonus relation, 1/£B% = £Bi which is relevant to the numerical 
evaluation of a range of integrals. A similar relation would map 3+1 integrals to 1+3. 
With the help of these relations we can reduce the propagator calculation to four converging 
integrals, B^(l), 6|(1, 1), <6|(2, 1) and 23^(2,1,1), the last of which vanishes for 3+1, and 
one infrared diverging integral, B^(2). 

For numerical calculation, setting n = rij and using the well-known Schwinger repre- 
sentation 



4l| . the integrals are rewritten as 
(_l)n roc 



B^(q, tlx, n 2 ,n 3 ,n 4 ) 



29-»r(g) 



dX \^e~ m2 ^ 2 J] 



d r ' 



dx r 



exp x I (x) 



(B9) 
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2 + 2 



3 + 1 



with 



1 0.15493339023106021408483720810745 0.15493339023106021408483720810745 

2 0.27309284159872576605248340464151 0.39746855267384293273601515359921 

3 0.34489265148380372504682107682350 0.63909144650207181551150685543306 

4 0.39476479072982726208009053119503 0.87599824804148552213829648851546 

5 0.43269826887073920594692758485023 1.10989753717092264667311285543088 

6 0.46325286965085016997908881189910 1.34197554748773397230022959566745 

7 0.48882308585244885903279712081054 1.57290464634371239724254250061935 

8 0.51080767194929084936132635888190 1.80307278894183022754245621178021 

9 0.53009088018361126653278475943989 2.03271449270324135251760022607330 

10 0.54726616387832592956739420406785 2.26197832054321458404301146488387 

TABLE VII: S ? (1) for 2 + 2 and 3 + 1 as a function of £ to 10~ 32 precision. 



a -exp^J (x) = ^r^(-ir fc M 3 ^4(x), (B10) 



dx 

k=0 

where M£' are partition multinomial coefficients and Ik(x) are modified Bessel functions of 
the first kind |47l l48j . Since we want to focus attention on the massless case (m 2 — > 0) and 
keeping only the non-vanishing terms either divergent or finite, power-counting shows that 
the integrals, B^(q, n 1; n 2 , n 3 , n<±) for q — n< 1, are infrared finite. As a result their value can 
be directly calculated from Eq. (|B9j) by setting m 2 = 0. Numerical evaluation to, say, 32 
digits becomes trivial with the use of a numerical integration package available with programs 
like Mathematica or Maple (see Table IVTij) . Furthermore, simple dimensional arguments give 
B^(q,nf,n c ) = 0((£log£) 9_ri/ /£) for q > rif and B^(q,nf,n c ) = 0{^ q ~ nf ) for q < rif for the 
2+2 case while for 3+1, B^(q,nf,n c ) = 0(^ q ~ nf ). From this, the asymptotic behaviour of 
the one-loop corrections can be easily obtained, by noting lim^oo jB|(2, l)/(^ 2 i3^(l)) = 1/6 
for 3+1. 

We can now described the computation of the relevant part of B^(q, ni,n 2 , n 3 , 0,4), q — n> 
2, which is infrared divergent. Using the asymptotic expansion for large x of I v {x) 

MJ V^^(2z)W(i/-Jfc + l/2)' 1 ] 

the leading and sub-leading behaviour is easily determined. This can be reexpressed in terms 
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of T(i, m 2 ), i = . . . q — n — 2 whose integral representation can be directly subtracted from 
the integrand, leaving an m 2 — ► converging integral. For example, in the case n — 0, the 
constants 6j are defined as the expansion of Iq(x) a , which are rational numbers multiplied 
by 7r -2 . The divergent part of B^(q) is then given by 

^ d ^ bi- 2 T(q-i) eVj 2 

; ; > — : Vt" log 771 . B12 

r(g) ^ 2 l (m 2 )9- 1 2PT(q) & V ; 

Following the literature, the finite contribution to 



B ( (2) = -f^(- logm 2 - 7 e + F (0) + 0(m 2 ) (B13) 



is calculated first and the finite part of the other integrals is defined up to it. Using 

r(0,m 2 ) = - \ogm 2 - lE -Y , (B14) 



nn\ 

n=l 



and taking the lim m 2^ , yields 

,2 r2 



+ jTdA(|;Ae X p-<«+^/nA)/;f(- 



For example 



^(3) = 3^ + tJ^ N™ 2 + 7^ - Fo(0) + /^(3), (B15) 

B|(3,l) = -|^(-logm 2 - 7s + F (0)+/^(3,l) (B16) 

i3 e c (3, 1) = ^ (-logm 2 - lE + F (0) + /£|(3, 1). (B17) 
which using Eq. IB5I and taking lim m 2^ m 2 £>g(3), must satisfy 

(4 - d)fB{(3, 1) + |/£|(3, l) + = 0. (B18) 

Numerical values for and Fq are presented for a range of anisotropies in Ta- 
bles f)VHIVHI[) . 
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2 + 2 



3 + 1 



1 4.3692252338747587180021767477 4.3692252338747587180021767478077 

2 3. 1818862274504847285875097585 3.0013807239614706751354255882375 

3 2.2336373613788360098511237787 2.0763302182432342962706829473432 

4 1.5519382904759206782491885601 1.4343462876751234663331108882378 

5 1.0353940883292892812839278453 0.9504439561174006940006958674003 

6 0.6235592513202110800048370509 0.5632462682606760634612041417967 

7 0.2823597812415587497126918746 0.2405570583083765791881447111856 

8 -0.0084789916018349072700020647 -0.036178339243301342993891362519 

9 -0.2617758427992475365501484308 -0.278533198254716148526058151986 

10 -0.4860763765593393395653061782 -0.494190516782956896533394268235 

TABLE VIII: F (g) for 2 + 2 and 3 + 1 as a function of £ to 10" 28 precision. 
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